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^^ ■ We use the quantum jump approach to study the entanglement dynamics of a 

\^ . quantum register, which is composed of two or three dipole-dipole coupled two- level 

^P \ atoms, interacting with a common environment. Our investigation of entanglement 

dynamics reflects that the environment has dual actions on the entanglement of the 
qubits in the model. While the environment destroys the entanglement induced by 
^ , the coherent dipole-dipole interactions, it can produce stable entanglement between 

d ' the qubits prepared initially in a separable state. The analysis shows that it is the 

^ i' entangled decoherence-free states contained as components in the initial state that 

contribute to the stable entanglement. Our study indicates how the environmental 
noise produces the entanglement and exposes the interplay of environmental noise 
^ ' and coherent interactions of qubits on the entanglement. 
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1 Introduction 



Entanglement is a remarkable feature of quantum mechanics and plays a fun- 
damental role in quantum computation and quantum information processing. 
Due to its fundamental importance for quantum information processing, it 
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has attracted much attention in recent years [T]. The entanglement can be 
produced either by direct interactions between qubits [2], or by indirect in- 
teractions between the qubits through a third party P]. However, both of the 
above processes are constrained to the closed system where the influence of 
the environment is neglected. In real situations, a quantum system can never 
be isolated and will inevitably interact with its environment. A severe effect of 
this unwanted interaction is decoherence which generally leads to the degrada- 
tion of quantum coherence and entanglement. Up to now, decoherence remains 
to be the main obstacle to the practical implement of quantum information 
processing. Recently, the influences of the environment on the entanglement 
have been investigated extensively [Hl5|6f7f8f9|10fllj . It was found that in cer- 
tain cases the environment can act as a third party and be constructive to the 
production of the entanglement, in contrast to its destructive role in the or- 
dinary situations. That is to say, depending on different coupling mechanisms 
between the environment and quantum system in different models, the envi- 
ronment has dual nature influenced on the entanglement between the qubits. 

In the present work, we will show that this dual influences can be present 
simultaneously even in one model under certain given coupling mechanism 
between the environment and system. For this purpose, we model a quan- 
tum register by two or three two-level atoms with dipole-dipole interactions. 
We study its entanglement dynamics (i.e. the time evolution of entanglement 
quantity [8]|12fl3fl4|[T5] ) under the influence of a common environment. The 
direct interactions between the qubits have not been considered in the previ- 
ous works [H[5]|5f7f5] . The environment used here is treated as a reservoir of 
quantum harmonic oscillators with inflnite degrees of freedom [16] , in contrast 
to the single-mode model used by Kim et al. [Ij. By introducing a collective 
mode consisting of all the qubit modes and eliminating the enormous irrele- 
vant degrees of freedom of the environment, a quantum master equation of 
the atomic system can be obtained under the Born-Markovian approximation. 
We use the quantum jump approach [TTj to study the time evolution behav- 
iors of bi- and tri-partite entanglement deflned by the concurrence [18] and 
the negativity [19] , respectively. Originally, the quantum jump approach [T7| 
was proposed to study the transition dynamics of a single atom system (for 
a review, see, e.g. [20]) and has been used to obtain the operator solution of 
the master equation [20PT] . The merit of using quantum jump approach is 
that it supplies a simple way to obtain the solution of master equation in the 
multipartite system with a very large Hilbert space. The dynamical behavior 
of the bi- and tri-partite entanglement shows clearly that the environment has 
dual nature influenced on the entanglement in the system studied. On the 
one hand, it destroys entanglement induced by the dipole-dipole interactions, 
which shows the decoherence nature of the environment. On the other hand, 
it induces stable entanglement in an incoherent way, which shows the third 
party role of the environment. Different from the previous works [H[5]|6ll71l8] . 
the dual nature of the environment is shown simultaneously in the same sys- 



tern in our situation. The stable entanglement induced by the environment is 
very different from the simple oscillation behavior of entanglement with time 
obtained by Kim et al. [4J, where the memory effect of the single-mode envi- 
ronment was considered. Here we point out that the environment in our model 
as a reservoir has no memory effect under the Markovian approximation. The 
environment induces entanglement in an incoherent way since with infinite 
degrees of freedom it introduces an irreversibility to the dynamics of the qubit 
system. 

The paper is organized as follows. In Section II we introduce the model and the 
quantum jump approach. In Section III we study the entanglement dynamics 
of a two-qubit system. The entanglement dynamics of a three-qubit system is 
discussed in Section IV. Finally, Section V is devoted to a brief summary. 



2 Model and quantum jump approach 



We consider a quantum register composed of A^ (=2 and 3 in the following 
discussions) identical two-level atoms coupled by dipole-dipole interactions. 
These atoms are assumed to be located very near, so that they feel a common 
environment and interact collectively with the environment. The environment 
is modeled by an electromagnetic field with infinite degrees of freedom which is 
assumed at vacuum state. If the collective interaction between atoms and the 
common environment is very weak and the characteristic time of correlation 
function of the environment is very short compared with the time scale of 
system evolution, then Born and Markovian approximations are valid. By the 
standard method of quantum optics tl6j, the master equation of the reduced 
density matrix of the system can be derived 

^ = -i[Hi, p (t)] + |(2J_p (t) J+ - J^J_p (t) - p (t) J^J_), (1) 
where 7 is the decay constant of the collective mode represented by 

N N 



J± = E^±^ J^=ll 



(y\ 



and 

N 

Hi = Y.^gij(r+(y- + h.c.) 

i>j 

describes the dipole-dipole interactions of the qubits. The first term on the 
right-hand side of Eq. ([1]) generates a coherent unitary evolution of the density 
matrix, while the second term represents the decoherence effect of the envi- 
ronment on the system and generates an incoherent dynamics of system. It 
should be noted that in general case the strengths of dipole-dipole interactions 



are dependent on the positions of the atoms. But in our model, the separation 
of the atoms are very small compared to the wavelengths of the most relevant 
field modes, so that the fields induce a global dipole-dipole interactions(in the 
long wavelength limit), i.e. gij = g. Thus the interaction Hamiltonian becomes 

/ A^ 

Hi = g[ J+J- - Jz- — 



2 
By using individual mode representation Eq. ([1]) can be rewritten as 

^ = -z[Hj, p (t)] + 1J: {2alp (t) ai - a^alp (t) - p (t) a^-). (2) 

It should be pointed out that the terms with i = j in the sum on the right- 
hand side of Eq. ([2]) denote the individual dissipation of each atom due to the 
environment, while the i ^ j terms describe the couplings between the atoms 
indirectly induced by the common environment. 

The master equation ([1]) has been used for many years to study the superra- 
diance effect involving the interaction of collective atomic systems with the 
radiation field [22|23] . The collective interaction between the atoms and the 
radiative field can induce two typical atomic states with an enhanced (the 
superradiant state) and a reduced (the subradiant state) spontaneous decay 
rate. When the atoms are confined in a very small region, the subradiant 
state is completely decoupled from the environment, and therefore can be re- 
garded as a decoherence-free state. A related experiment on the superradiant 
spontaneous emission of two trapped ions has been reported [2] 



The common feature of the quantum master equations is the existence of the 
sandwich terms, where the reduced density matrix of the system is in between 
some quantum excitation and de-excitation operators. In general, the master 
equation is converted into c-number equation [TOf^ . In this paper, we use the 
quantum jump approach to obtain the exact operator solution of the master 
equation ([1]). As mentioned above, the quantum-jump approach was proposed 
to describe single experimental realizations of simple quantum systems (i.e. 
two- and three-level atoms) [17] and has been extensively used in quantum 
optics [20]. By using the quantum jump approach, it is straightforward to 
obtain the exact operator solution of the master equation. The advantage of 
using this method is that one does not need to deal with each element of 
the reduced density matrix individually, the time-dependent solution of the 
reduced density matrix can be obtained directly from the formal solution under 
the actions of operators. 

By the quantum jump approach the reduced density matrix p (t) can be writ- 



ten as a sum of the conditional density matrices p^"-\t){n = 1, 2, ■ ■ ■ ) as follows 
Pit) = Y.P^'Ht) = T.TrF{QnPtotm, 



n 



where ptot is the total density matrix of the system and the quantized envi- 
ronment, TrF{-} is the partial trace over the modes of the environment, and 
Qn is the projection operator on the state of the quantized environment field 
that contains n photons. The p*^") (t) denotes the reduced density matrix of 
the system with n photons detected in the environment. The time evolution 
equation of p*-"^ (t) can be derived readily, 

^^^ = -^{H^ffP^-'> (t) - pW (t) ny + ^J_p(-'^ (t) Ml - <5„o), (3) 

where Heff = -ff/ — yJ+J-. Eq. (3) reflects that the time evolution of the 
reduced density matrix consists of smooth evolutions which are interrupted 
by instantaneous quantum jumps. It is clear to see that the hierarchy of the 
conditional density matrices p*-"-'(t) is terminated by instantaneous quantum 
jumps automatically. Due to the dissipative effect of the environment, the 
effective Hamiltonian -f^e// governed the smooth evolution is a non-Hermitian 
operator. It is straightforward to write down the formal operator solution of 
Eq. m 



pW(t) =7 rt/t'e-^^^//(*-*')j_p("'i)(t')J+e</(*-*'\ (4) 



where p(0) is the initial state. Thus, the total reduced density matrix of the 
system can be obtained by the summation over all the conditional density 
operators p*^"-* (t). In the following we explicitly consider the two-qubit and 
three-qubit systems. 



3 Entanglement dynamics in a two-qubit system 



To simplify the following calculations, we use the collective state representa- 
tion [23j. The representation space is spanned by the common eigenstates of 
the complete set of commuting operators {J^, Jz} for two-qubit system. This 
collective state representation can also be obtained by the angular momentum 
additive rules and the Clebsch-Gordan coefficients from the product state rep- 
resentation of the individual ones |25j. The two-qubit basis denoted by \j,m) 
includes the spin singlet and triplet 




Fig. 1. The time evolution of the concurrence C{t). The pair of atoms is initially 
prepared in the state p{0) = \ -\ — )(H — | for (7 = 1.0 (solid line), 0.2 (dashed line), 
and (dotted line) and 7 = 0.1. 

|o,o) = -^(| + -)-!-+)) 
|i,o) = -^(| + -) + !-+)), 

|1,-1) = |--), (5) 

where |±) are the eigenstates of o"^ with elgenvalus ±1, respectively. The J± 
act on the collective basis by 



J±\j,m) = ^j{j + 1) - m{m ± l)|j, mil). 



When a total system is composed of two subsystems described by a two- 
dimensional Hilbert state, the general measure of entanglement between these 
two subsystems is the entanglement of formation [26]. This quantity can be 
analytically calculated as a function of the concurrence C [18], which can also 
be taken as a measure of the entanglement. The concurrence is defined as 



C = max(0, Ai — A2 — A3 — A4), 



(6) 



where Aj(i = 1, ■ ■ ■ ,4) are eigenvalues of the matrix [p^^'^{cry (g) (Ty)p*{ay ® 
^y)P^ Y ■ Here p* is the complex conjugation of p. The concurrence varies 
from C = for a disentangled state to C = 1 for a maximally entangled state. 
We shall use the concurrence to quantify the degree of entanglement in the 
two-qubit system. 

In the following numerical investigations, our guide line is to explore the fol- 
lowing topics of entanglement dynamics: 



(1) Two kinds of entanglement generation mechanisms and their different roles 
played: direct coherent generation due to the atomic dipole-dipole interaction 
and indirect incoherent generation due to the common environment coupling 
to the qubits. From Eq. ([T]) we can see that the atomic dipole-dipole interac- 
tions in Hi produce a dynamical coherence which oscillates in time, and the 
common environment induces an effective interaction which produces the en- 
tanglement in steady states. Besides, the common environment also produces 
a dissipative term which makes all the eigen states of the operator of the mas- 
ter equation decay and finally die out except the steady states with zero eigen 
values [H]. Thus in the dissipative dynamics, the dynamical coherence as a 
transit phenomenon will be destroyed finally by the dissipation, only the envi- 
ronment induced coherence and entanglement in the decoherence-free steady 
states survive. In the following numerical study, the two kinds of entanglement 
generation show clearly in the figures. 

(2) How does the environment produce an entangled steady state from a sep- 
arable initial state? This issue is closely related to the generation of an entan- 
glement source and is thus very significant. 

(3) How does the entanglement measure evolve and change in time from an 
initial entangled state under the influence of the common environment? This 
issue is related to the decoherence of a multi-qubit system and is thus practi- 
cally helpful in the understanding of the decoherence dynamics. 

(4) As the decoherence-free subspace has more than one dimensions, the steady 
state will be a probabilistic mixture of the decoherence-free states. How does 
one choose a proper initial state to get a desired steady state? 

In order to show explicitly the entanglement dynamics of the two-qubit system, 
we take the initial state p(0) = | H — )(H — | as our first example. Expanding 
this state in the collective basis of Eqs. ([5]) and substituting it into Eqs. (jl]), 
we can obtain analytically the solution of the master equation ([T]) after some 
calculations 

g-27t I _ g-27t 

pW = ^-|1,0)(1,0| + ^ |i,_i)(i,_i| 

1 p—i'igi—'yt pi'igi—'yt 

+ -|0,0)(0,0| + ^ 11,0)(0,0| + ^— 10,0)(1,0|. (7) 

The time evolution behavior of the entanglement measure, i.e. the concurrence, 
can be determined by 

C{t) = iy^(e-27t - 1)2 + 4e-27t(sin2c/t)2. (8) 

From Eqs. (^^ and (^ we can see that the steady solution of the atomic system 




Fig. 2. The time evolution of AC = C\ 
7 = (dashed hne). 



9=1.0 



C|g=o with 7 = 0.1 (sohd line) and 



under the long time limit is 



Ps 
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-1)(1,-1| + -|0,0)(0,0| 



which is an equal probabilistic mixture of the extremal states |j, —j) and an en- 
tangled mixed state with Cg = 0.5. The extremal states \j, —j) is decoherence- 
free states, as shown in Ref. [27]. One of the decoherence-free states, i.e. |0,0) 
is the so-called subradiant state [22] and maximal entangled pure state. This 
means that due to the symmetry of the interaction between the system and 
the environment, the system induces a collective entangled state which is im- 
mune to the decoherence. This entangled decoherence-free state contributes to 
the stable entanglement produced by the environment. From the purification 
scheme proposed in Ref. [26j many pairs of this states can be used to distill 
out a maximal entangled state |0,0) with probability 1/16. 

The time evolution behavior of the concurrence is showed in Fig. [H We con- 
sider three cases of the dipole-dipole interaction: absence {g = 0, dotted line), 
weak {g = 0.2, dashed line), and strong {g = 1.0, solid line). In the case of 
g = 0, the entanglement increases monotonously and approaches to a stable 
value. Since there are no dipole-dipole interaction between the qubits, it is 
quite clear that this entanglement is induced purely by the environment. This 
is a typical example that the environment plays a constructive role in the 
entanglement formation between qubits. When the dipole-dipole interaction 
is switched on, the entanglement induced by the environment is not affected 
but an oscillating dynamical contribution added to the total entanglement is 
observed. The stronger the dipole-dipole interaction is, the faster the entan- 
glement oscillates, as seen easily from Eq. ([8]). To see the contribution of the 
dipole-dipole interaction to the entanglement, in Fig. [2]we plot the quantity of 
AC = C\g=i,o — C\g=o as a function of time. It is shown that the entanglement 



induced by the dipole-dipole interaction decays rapidly as shown by sohd hne, 
which indicates that the environment plays a destructive role to the entangle- 
ment induced by the coherent interaction. To confirm this, we switch off the 
environment, i.e., 7 = 0. In this situation, the entanglement induced by the 
coherent interaction oscillates without dissipation, as the dashed line shown 
in Fig. [21 To summarize, the environment can act as a third party to induce 
the entanglement between the qubits incoherently. It also plays an usual role 
of dissipation to destroy the entanglement induced coherently by the dipole- 
dipole interaction. The dual nature of both construction and destruction roles 
played by the environment on the entanglement generation can be seen clearly 
by the above numerical analysis. 

To explore the dependence of the steady state entanglement on initial states, 
in the following we take the Werner states as the initial states. The Werner 
states read 






where I4 is a 4 x 4 identity matrix, |\E'+) = |1, 0) and |^-) = |0, 0) defined by 
Eq. ([5]), and p G [—1/3, 1] denotes the fidelity of W± to \'^±)- It is known that 
W± are entangled states for p > 1/3 and C± = ^^-^. In a similar way, one can 
obtain the time evolution behavior of p± as 



P+{t) = ^ 11,1)(1,1| + { ^ + ^ > |1,0)(1,0| 

^|0,0)(0,0|+ ^ ^'^ ^—^ — ^ ' ^^ |1,-1)(1, 



^(l^p)t^:^ni±7Sl|i,_i)(i._i|,i±^|0,0)(0.0|. (10) 

Because only the diagonal forms of collective states are involved in the Werner 
states iQ, the coherent interaction Hj of Eq. ([1]) has no actions on the Werner 
states. 

The concurrence can be calculated in a similar way. Before discussing the 
dynamical behavior of the concurrence, it is useful to consider the long time 
behavior of the concurrence. With t — > 00, p±(t) become 




Fig. 3. The time evolution of the concurrence C+ with different p. The pair of atoms 
is initiahy prepared in the state p{0) = W^ and 7 = 0.1. 

P_. = i±*!|*.)(*_H-i^|l.-l)(l.-l|, 

with the steady entanglements as C+s = ^-^ and C^g = "^^^5 respectively. 
The above results show that different initial states lead to different steady 
states with different entanglements. Also from the limit case with p = 1 one 
can see that \l/+ is sensitive to the decoherence, while \E'_ is not. 

Below we discuss entanglement dynamics in the general case. Fig. [3] and Fig. H] 
show the concurrences of p±(t) as functions of p and t. When p < 1/3, a com- 
mon characteristic of p±{t) is that they have no entanglement initially. With 
time a stable entanglement is formed through the environment. When p > 1/3, 
the concurrences of p±(t) show apparently different time evolution behaviors. 
For p+(t), the entanglement contained in the initial state W+ is firstly de- 
stroyed by the environment due to the decoherence effect of the environment. 
The entanglement experiences a "sudden death" [T3]. After a finite time du- 
ration the entanglement is reconstructed due to the constructive effect of the 
environment. Physically, the initial contribution to the entanglement in W+ 
comes from the component of |^+). With time, |\1'+) (i.e. |1,0)) approaches 
to the product state |1, —1) and its entanglement disappears gradually. In this 
process, the decoherence- free state |\E'_) is unchanged and the corresponding 
entanglement gradually dominates the entanglement of p+{t) and finally, it 
reaches the steady state value in the long time limit, i.e., (1 —p)/4, as men- 
tioned above. Different from p+(t), the initial entanglement in p-(t) does not 
decay. This is because the main contribution to the initial entanglement in this 
case comes from the decoherence-free state I^E'-). The entanglement induced 
by the environment is superposed on the initial entanglement. When p = 1, 
the Werner state reduces to the pure state I^E'-) which is independent of time 
and free from decoherence. 
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Fig. 4. The time evolution of the concurrence C_ with different p. The pair of atoms 
is initiahy prepared in the state /7(0) = W- and 7 = 0.1. 



From the special and general analysis above, one sees that the entanglement 
of the steady state induced by the environment is rooted in the seeds of the 
entangled decoherence-free states which as components have already been con- 
tained in the initial states. This property can be used in choosing a proper 
initial state to obtain a desired entangled steady state. 



4 Entanglement dynamics in a three-qubit system 



The three-qubit basis can be constructed from the addition of angular mo- 
menta Ji, J2 and J3. The representation space is spanned by the common 
eigenstates of the complete set of commuting operators {JI21 J^i Jz}- The col- 
lective basis |ji,2, j, TT^) can be obtained by virtue of the corresponding Clebsch- 
Gordan coefficients 
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1.2.2> = ! + ++>■ 

ii-?) = l--) 

'2' 2' ' '■ 

1, 2- 2> = 7j(l + +-> + I + -+> + I - ++». 

The J± act on the collective basis as follows 

J±\ji,2d,m) = ^j{j + l) -m(m±l)|ji,2,j,m±l), 

which are independent of the quantum numbers ji^2- 

To discuss the entanglement dynamics in the three-qubit system, we use the 
negativity proposed by Vidal and Werner to quantify the degree of entan- 
glement |l9j. The idea of this measure of the entanglement comes from the 
Peres-Horodecki criterion for the separability of bipartite systems [28]. The 
negativity was originally introduced to an arbitrary two-qubit state p and 
defined as [T91I29] 

iV(p) = -2^^r 
■/ 

where /i~ is the negative eigenvalue of the partial transpose of p with respect 
to the i—th qubit, i.e. p^\ Given a bipartite state one can calculate the partial 
transpose p^* of the density operator. The state is exactly separable if p-^' is 
again a positive operator. However, if one of the eigenvalues of p^* is negative, 
the state is entangled |28]J. From this viewpoint, the negativity is used to 
quantify the degree that p^' fails to be positive and to represent the strength 
of quantum correlation between the two subsystems. It was proved that the 



relation between C{p) and N{p) is N{p) ^ ^(1 - C{p)y + C^{p) - [1 - C{p)] 
for any two-qubit state p9] . The merit of using the negativity to quantify the 
entanglement is that it allows us to investigate the entanglement properties 
between part i and the sum of other components in the multipartite system. 

Similarly, we choose an explicit initial state, e.g., p(0) = | H \-){-\ 1-| to 

discuss the entanglement dynamics in this case. Expanding this initial state 
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Fig. 5. The time evolution of the entanglement measured by the negativity Na when 

the three-atomic state is initially prepared in the state p(0) = | H |-)(H 1-| for 

g = 1.0 (solid line), 0.2 (dashed line), and (dotted line). The decay constant is 
chosen to be 7 = 0.1. 



in terms of the coupled basis of Eqs. flTTl) 



p(o) = lx)(xl, 

'^5/ 2 ' 2 2' 2' 2 2' 



:i2) 



and using the quantum jump approach, the time-dependent solution of the 
master equation can be obtained analytically 
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In Fig. [5] and Fig. [6] we show the time evolution behavior of the negativity 
A^^ and Nb corresponding to the partial transpose with respect to the first 
atom A and second atom B, respectively. In the present case, Na = Nc, 
where C denotes the third atom. One can see that the system shows a similar 
property of entanglement production as the two-qubit system. In absence of 
the direct interactions, the entanglement induced by the environment increases 
monotonously, then approaches to a stable value at a large time scale for both 
Na and Nb- After the interactions are switched on, the oscillation occurs. 
Moreover, during the evolution of time, this oscillation is suppressed gradually 
by the environment. From Eq. flT^ it is readily seen that the time-dependent 
solution asymptotically tends to the steady state 



M \/i I 1/ \// I 

P.-jll.j.^Xl.j.l-l + gW'^X^I 

which is a probabilistic mixture of decoherence-free state |1, |, ^) and |^^) 
[27] . In the present three-qubit case, there are two entangled decoherence-free 
states which contribute to the entanglement of the steady state, i.e. |0, |, ^) 
and |1, |, ^). It is interesting to note that pg is not a simple mixture of three 
decoherence-free states. It still involves quantum coherence between |0, |, ^) 
and |1,|,^). By comparison of the initial state Eqs. (fT2l) with the steady 
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Fig. 6. The time evolution of the entanglement measured by the negativity Nb when 

the three-atomic state is initially prepared in the state p(0) = | H \-){-\ 1-| for 

g = 1.0 (solid line), 0.2 (dashed line), and (dotted line). The decay constant is 
chosen to be 7 = 0.1. 



state Eq. (fl^ . we see that although the jump has occurred, i.e. from \ipi) 
jumps to 1^^) by the action of J_, the superposition coefficients between 

|0, ^,jz) and |1, ^,jz) are not changed. This is because that the subspaces of 
|0, |, jz) and |1, |, J2) are degenerate under the actions of both the coherent 
interaction Hj and the dissipative operators J±. Anyway, the steady state ps is 
still an entangled state due to the entangled decoherence-free states contained 
in \ip^), which is a linear superposition of the two entangled decoherence-free 
states. 



The above discussion in the three-qubit case further confirms the conclusion 
that the environment has dual nature on the entanglement. Different from the 
two-qubit case, the Hilbert space of the three-qubit system reduces into three 
decoupled subspaces |1,|,J2), |0, |,jz), and \l,^,jz), each of which has its 
own decoherence-free state [27]. Only two of the decoherence-free states are 
entangled. Under the time evolution, each component in the three subspaces 
of the initial state approaches its corresponding decoherence-free state. The 
necessary condition to induce stable entanglement is that the initial state 
contains the components in the two subspaces with entangled decoherence- 
free states, i.e. in the subspaces of |0, hjz) and |1, ^,jz)- 



5 Summary 



In conclusion, we have investigated a model of quantum register, which is com- 
posed of two or three atoms and coupled to a common environment. Using the 
quantum jump approach, the time-dependent solution of the master equation 
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and the entanglement dynamics of the system are studied analytically and nu- 
merically. The dual nature of the common environment on the entanglement 
and its dissipative dynamical origin are explored in detail based on the eigen 
solutions of the operator of master equation. On the one hand, due to its in- 
duced dissipative term in the master equation, the environment destroys the 
entanglement induced coherently by the atomic dipole-dipole interactions. On 
the other hand, due to its induced effective atomic interactions in the master 
equation, the environment can incoherently induce the entanglement among 
qubits in the decoherence-free space. The entanglement dynamics studied in 
the present paper addresses the constructive role of the environment on the 
entanglement production and suggests to make use of the above positive role 
to construct an environment-assisted entanglement production in quantum 
information processing. 
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